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Proton decay is a generic prediction of GUT models and is therefore an important channel to detect 
the existence of unification or to set hmits on GUT models. Current bounds on the proton lifetime 
are around 10^^ years, which sets stringent limits on the GUT scale. These limits are obtained 
under 'reasonable' assumptions about the size of the hadronic matrix elements. In this paper we 
present a non-perturbative calculation of the hadronic matrix elements within the chiral bag model 
of the proton. We argue that there is an exponential suppression of the matrix elements, due to non- 
perturbative QCD, that stifles proton decay by orders of magnitude - potentially 0(10"^"). This 
suppression is present for small quark masses and is due to the chiral symmetry breaking of QCD. 
Such a suppression has clear implications for GUT models and could resuscitate several scenarios. 



I. INTRODUCTION 

Proton decay is an important prediction of numer- 
ous models of beyond the standard model (SM) physics. 
The decay is induced by higher-dimensional, baryon- 
number violating operators, suppressed by some high 
scale, Mgut Although there are considerable un- 
certainties in the strength of the baryon violating cou- 
plings, they are usually assumed to be 0(1). One can 
then make an order-of-magnitude estimate of the pro- 
ton lifetime. For example, a, p ^ n + £ decay channel 
contributes 



2^|A|2|(^|0|p)|VK 



(1) 



where |Ap contains the non-hadronic part of the ma- 
trix element along with any perturbative coefficients, e.g. 
A ~ 1 /Mqjjj. when baryon- number violation is caused by 
four-fermion interactions. For a more complete estimate 
of Eq. 0, one also needs to know the hadronic matrix el- 
ements, which have previously been estimated in lattice 
calculations [2H5i chiral lagrangians [5] and other ap- 
proaches. In this paper, we calculate the hadronic matrix 
elements within the chiral bag model of the proton [5HT5] . 
The chiral bag model is the marriage of two interesting 
phenomenological descriptions of protons. At short dis- 
tances the proton is described as a 'bag' of free fermions 
- massless, or nearly massless fermions in a spherical 
bag subject to physically motivated boundary conditions. 
Outside of the bag radius, the proton is described with a 
Skyrme [HIH] model. 

Combining the two descriptions provides an improved 
model of the proton. Unlike the original bag models 
(the so-called MIT bag [16U21j). many physical results 
in the chiral bag description are insensitive to the bag 
radius P^24j . Simply, the smaller (larger) the bag, the 
more the proton is carried in the Skyrmion (bag). If 
the bag radius is taken completely to zero, we recover a 
purely Skyrmionic description of the proton. The Skyrme 
model has had considerable success in describing some 
baryonic properties |14j il5.i (25, 2^ , however it is inade- 
quate for describing proton decay since baryon number 
is identified with the winding of the pion field into the 



Skyrmion and is therefore topologically conserved. How- 
ever, for a Skyrmion with a 'hole', topological conserva- 
tion is not exact and the configuration can be un- wound, 
a process which can be interpreted as the Skyrmion (pro- 
ton) decaying to the topologically trivial ground state. 
The intuitive expectation is that the unwinding of a 
Skyrmion involves tunneling through a potential barrier 
and therefore comes with some exponential suppression. 
The smaller the hole, the harder it is for the proton to 
shed its topological portion and decay, and as a result 
there could be some topological suppression of the near- 
Skyrmion-to-vacuum transition. The chiral bag model is 
exactly a Skyrmion with a hole, so it is an interesting 
laboratory to study proton decay. Does the intuitive pic- 
ture or proton decay in the chiral bag model hold up? 
How large of a suppression is there? How sensitive is it 
to the bag radius? These are the sort of questions we 
aim to address here. 

The setup of this paper is the following: we introduce 
the basic facts of the chiral langrangian and the Skyrme 
followed by an introduction to quark 



nil. Section IV contains a simplified 



solution in Sec. [H 
bag models (Sec. 
calculation of Skyrmion unwinding where we neglect the 
interior (bag) dynamics. This calculation is subsequently 
improved in Sec. |V] with a detailed calculation of the 
bag energy in the presence of time-dependent boundary 
conditions. Numerical results are presented in Sec. |IV| 
and Sec. |Vl We end with a discussion of our results and 
some directions for future work. 



II. CHIRAL LAGRANGIAN AND SKYRME 
SOLUTION 

The low-energy effective theory for QCD is given by 
a non-linear sigma model, describing pions as the Gold- 
stone bosons from the spontaneous symmetry breaking of 
chiral symmetry. In our setup, we work with two mass- 
less quark flavors, giving rise to an exact chiral SU{2) 
symmetry. The pion Lagrangian is 



= ^Tr [d^U^d'^U] 



(2) 
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where U is the SU{2) valued pion field. This is a four- 
dimensional non-linear sigma model, where the ground 
state is given by constant field that we are free to choose 
to be [/ = 1. 

It is convenient to introduce — Ud^Jj'^ = —Xj^^ 
which is an element of the su(2) Lie algebra. The total 
energy of a static field configuration is given by 



Eo = ^ I d^x Tr 



(3) 



Finite energy solutions must have J7 — 1 at infin- 
ity, hence finite energy solution can be compactificd to 
maps of S'^ — > SU{2) = S^. These maps are clas- 
sified by TTs^S^) = Z, so there exist topological non- 
trivial field configurations. However, via simple scaling 
arguments, Derrick's theorem [37] shows that these non- 
trivial field configurations are unstable against scaling. 
Skyrme [TH [TS] added an extra term to the Lagrangian 



(4) 



which is the unique, lowest dimensional, higher-order 
term that satisfies all the symmetry constraints and is 
second order in time derivatives. This term stabilizes the 
field against scaling, thus allowing for stable non-trivial 
solitons. 

For our purposes, it is convenient to Wick rotate and go 
to dimensionless variables. Specifically, we set x^^ ^ Ry^ 
where we have introduced the characteristic length scale 
R = , so from now on everything will be dimension- 
less. With this, the (Euclidean) Skyrme action becomes 



1 



d^y -TviXaXa] 



TY[Xa,Xb]^ 



(5) 



The symmetry group of this model is SU{2) l x SU {2)]^ 
where U — )■ LUR^, in perfect correspondence with the 
QCD symmetry group. The axial U{1) symmetry is 
broken because U is restricted to SU(2), which corre- 
sponds to the anomalous breaking of U{1)a in QCD. 
However, the baryonic U{l)v symmetry is not present, 
as J7 — e'"^C/e~'"^ = U, so there is no Noether current 
associated with the baryon current. However this theory 
has an extra conserved current, 



247r2 



Tr [XjyXaXfj 



(6) 



This current is conserved identically without invoking the 
equation of motion. This current is purely topological 
and its conserved charge 



B 



d^xB^{x) 



(7) 



counts the topological winding of the U field. Identifying 
B with baryon number, the solitonic solutions of the the- 
ory are interpreted as baryons. With this identification. 



baryon number is a topologically conserved quantity, so 
protons cannot decay. This property, and the extra stabi- 
lizing feature which it requires, makes it reasonable that 
proton decay may be suppressed more than naively ex- 
pected, even outside the strict confines of the Skyrme 
model. 

The proton is the stable solitonic solution (Skyrmion) 
with winding number one. It can be found by making 

the ansatz 



U — exp [i F{r) x ■ a] 



(8) 



supplemented by the boundary conditions that F{oo) — 
and F{0) = w. With these boundary conditions, U 
wraps around the SU{2) exactly once. For this ansatz, 
the energy is given by 



E\F] 



vr 
'2^ 



dr r 



,o sin2(F)\ 
F'^ + 2 



4 sin (F) /sin2(F) 



2F'- 



(9) 



Minimizing E[F] with respect to the function F deter- 
mines the profile of the Skyrmion, which is shown in 
Fig.[l] 
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FIG. 1: Profile of the Skyrmion solution. 

The boundary condition at r — > is set by the re- 
quirement that the Skyrmion has baryon number equal 
to one. Integrating the topological charge for the Skyrme 
solution from tq to infinity, we get 



B=-{F{ro)-lsm2F{ro)). 



(10) 



For F{0) — IT, the baryon number B = B{0) = 1 . 

Beautiful as this model is, experiments have clearly 
shown that the quark picture is the correct picture for 
small distance scales. The chiral bag model incorporates 
this by punching out a hole in the Skyrmion and replacing 
physics in that hole with a bag of free quarks; in essence, 
the hybrid bag models factorize QCD in short- and long- 
distance regimes. We will return later to how the hole 
effects Skyrmion properties. However, first we will review 
some properties of quark bag models. 
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III. CHIRAL BAG 

As a first approximation, the inside of the bag sim- 
ply contains free quarks. As such, it is described by the 
Minkowski path-integral 



DipDip exp • 



(11) 



where integration is restricted to a spherical region of 
radius rf,Qg. The boundary conditions, at r^ag, are chosen 
such that they respect the symmetries of QCD: 



{upr + u^PL)i' = (JsV- 



(12) 



Here is a radial unit vector pointing outwards, and 
the matrices U are the non-linear pion field at the bag 
boundary. One can now easily see that this model has 
the right SU{2)l x SU{2)r symmetry structure, of ijjR 
Ri'R, V'L LipL, U — >■ LUB) . For now we take U = 
exp[i0n • cr], = exp[z6'7^n • u], for some value of 9. 
However eventually we will make the identification 6 = 
^('"bag), the solution to Eq. (j9|. This boundary condition 
will then link the properties inside the bag to the Skryme 
solution on the exterior. The conditions in Eq. ( 12 ) imply 



no vector currents flow through the boundary, as 



thus 



(13) 



(14) 



Therefore, baryon and color currents through the bag 
boundary are zero"'^. However, unlike the original (MIT) 
bag models, the axial SU{2) current is continuous along 
the boundary of chiral bags once we identify 6 = F{r), 
because of the full presence of the axial symmetry. 

We have shown that the baryon number is confined for 
fixed 9. If 9 varies, the baryon number of the Skyrmion 
changes, but so does the baryon number of the bag! The 
Dirac sea eigenstates energies are modified by and depend 
nontrivially on the boundary rotation 9. Through the 
n ■ a term on the boundary, spin and isospin are linked 
together, so eigenstates must be classified by the sum 
S + I, rather than spin or isospin individually. One can 
show using the map ipE — >■ 75'0£ that the spectrum is 
invariant under E E and -> 6* + tt hence periodic 
with period tt. Using ipE ^ once can show that 



^ Because the boundary couples flavor states but leaves color un- 
touched and we neglect full QCD inside, there is a TVc-fold de- 
generacy of all eigenstates. Therefore we focus on one color eigen- 
state and multiply by Nc where appropriate in the rest of the 
paper. 



the spectrum obeys the symmetry E —E and 9 ^• 
TT — 9. Therefore the spectrum has a zero-mode at 9 — 
tt/2. The ^-dependent eigenvalues have been worked out 
in Ref. [TT| and sketched in Fig. |2] However, we are 
interested in properties of the bag as a whole, such as its 
total energy density or baryon number, rather than the 
individual eigenvalues. 




FIG. 2: Visualization of the eigenstates as function of 6. The 
spectrum at ^ = tt is the same as 9 — 0, but all levels go down 
by one step, n — n — 1. At 6' = tt/2 there is a zero-mode. 

The baryon number of the vacuum is defined with re- 
spect to a state in which every state is 'half filled. Every 
empty state therefore counts as baryon number — | , and 
every filled state as ^. Summing over all modes. 



i^sgn(£;„ 



(15) 



which must be suitably regulated. In Ref. [9], the sum 
in Eq. ( 15 ) was evaluated, and the baryon number for a 
chiral bag of radius r was shown to be 



1 



9 < 7t/2 
9 > tt/2 



1 



sin(20(r)) . (16) 



Identifying 9 with the value of the Skyrme profile func- 
tion at the boundary F(rhag) and summing Eqs. (16) 
and (10 1, aU dependence on rtag drops out and the net 
(Skyrmion -|- bag) baryon number remains constant at 
one.^ This shows that, indeed, the topological current is 



^ Except when 9 becomes less than 7r/2. We comment on this in 
section ^ 
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the baryon current and baryon number is a globally con- 
served current. This remarkable property is referred to 
as the Cheshire Cat principle [55] in the literature and it 
gives confidence to the idea that (at least some properties 
of) baryons can be consistently factorized into short dis- 
tance quark (and gluon) physics and long distance pion 
physics. In addition, the full SU{2)]^ x SU{2)ji symme- 
try is maintained, making the chiral bag model a model 
for baryons. 

Having seen that the quark bag and Skyrmion scenar- 
ios can be stitched together to form a consistent picture 
of the proton, we now want to study proton decay within 
this setup. We will proceed in steps. We first calcu- 
late the unwinding of a 'punctured' Skyrmion neglecting 
changes to the bag interior. Then we add two crucial 
ingredients, the Casimir energy of the interior bag and 
source terms for baryon-number violation in the action. 



Plotting the potential V^(A), there is a meta-stable mini- 
mum at A = 1, and an absolute minimum at A — )■ oo. Be- 
tween the minima is the potential barrier associated with 
unwinding the Skyrmion. Transmission through the bar- 
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IV. NEGLECTING THE INTERIOR 

As stated earlier, a punctured Skyrmion is not a topo- 
logically conserved state, so it can decay into the vacuum 
state. Our approach to the classical solution describing 
the decay is the following: we take a Skyrmion with fixed 
puncture size Vfjag, then consider dilations of that solu- 
tion, F{r) — T' F{X{t)r). The dilation parameter depends 
on time, with A(— oo) = 1. As A grows, more of the 
Skyrmion is sucked into the hole, and as A — > oo the 
Skyrmion disappears completely and we are left with the 
trivial ground state. Changing variables in Eq. ([5| to 
r — A(t)r, and plugging in the Skyrme ansatz, we find 



2A 



2 VA 
+ A 



A3 



'2sin^(F)+.^F- , .f^^!AFl+SsinHF)F- 



(17) 



where primes indicate derivatives with respect to r and 
dots indicate a time derivative^. Substituting the Skyrme 
solution F{r, t) into the above and performing the spatial 
integration, we are left with a one-dimensional problem. 



dT X"^ K{X) + V{X). 



(18) 



Because the integral for K{X) is convergent when one sets 
Tbag — in the lower limit, we have approximated K(X) 
by 



K(X) 



A B 
A3 ^ A5' 



(19) 



FIG. 3: The potential ^(A) for a Skyrmion with bag radius 
T'bag ~ 0.1 R, where we have subtracted off V^(l) = Vq. The 
tunneling probability is calculated by inverting the potential, 
V ^ —V and determining the bounce solution. For this par- 
ticular example, the turn-around point is A* ~ 28. 

rier can be calculated using instanton/WKB techniques, 
by which we construct a 'bounce' solution that starts at 
A(— oo) = 1, rolls down the inverted potential to the other 
side reverses, and ends up at A(-l-oo) = 1 (in Eq. ( [T7| , 
( jTs] ) we have already been working with the Euclidean 
action). The action of this 'bounce' solution is** 



5, 



tunnel 



A* 



dX ^KiX){ViX)-V{l)) 



(20) 



where A* is the turnaround point for the bounce solution. 
The width of the decay of the punctured Skyrmion is 
proportional to e~'^'"""='. 

The 'bounce' solution can be thought of as the 
Skyrmion unwinding, then rewinding. To get an idea for 
the size of the suppression, we can plug in some numbers. 
For puncture sizes of r{,ag/R — 0.05, 0.1, 0.3 we find expo- 
nential suppression of 5.7 x 10^^, 1.1 x 10^"*, 2.5 x lO^'^ 
respectively. From the height and width of the potential 
in Fig. |3j one may have expected a larger suppression. 
The reason the suppression is not larger is because of the 
A~3, A~^ terms in K{X); these terms quickly shrink as A 
increases, resulting in a smaller than expected tunneling 
action. 

The calculation at this point clearly depends strongly 
on the size of the puncture, an issue we will return to 
soon. However, by neglecting the bag interior, this cal- 
culation is missing some physics. As the Skyrmion is 



As we have gone to dimensionless variables, time derivatives are 
taken with respect to t = t/ R. 



* Note there has to be at local minimum at A = 1, unlike in Figjs] 
However just imagine a small change to make it a local minimum. 
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pushed into the hole, the value 9 = F{rbag A(t)) of the 
Skyrmion at the bag boundary changes. The changing 
boundary condition shifts the energy levels of the bag 
fermions, an effect we need to incorporate. 



V. INCLUDING THE INTERIOR 

To include the effects of the fermions in the bag we 
need to consider the effects of the fermionic path-integral. 
The resulting functional determinant, for static bound- 
ary conditions, was derived in Ref. [21] and is also 
presented in Appendix |A] Schematically, the fermionic 
path integral results in a ^-dependent functional de- 
terminant which, once regulated, can be shown to be 
exp[— T£'cas(^)], where Ecas is the Casimir energy. Since 
it depends on 9, the interior Casimir energy acts as an 
additional potential term in the action for 9, 



The Casimir energy is 



D[72?V'^2?Vexp I - y d^x{LE{9) -f i'^dr + H{9))^Ij^ } 
JdUcxp^-J d\LE{9) + Nc jdTE,as{9)^, 



(21) 



where the Hamiltonian is H{9) ~ —17"'^ together with 
the boundary conditions (Eq. [12] with angle 9). We have 
reverted back to using 9 to describe the boundary angle 
in this section, however the reader should keep in mind 
that 9 is set by the Skyrmion solution at the interface, 
9 = F{riyag)- The Casimir energy of the spectrum of 
the Hamiltonian is defined as Ecas — —\^n \^ri\- We 
inserted a factor of Nc to account for the color degener- 
acy. The total energy of a static solution is thus given 
by the sum of Skyrme energy (Eq. ([9])) and the Casimir 
energy (Eq. ^). In Fig. [Ij the total energy is shown 



and one can see that it is remarkably flat as a function of 
the bag size, further evidence of Cheshire Cat principle. 
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FIG. 4: The energy of the Skyrmion (dashed, in blue), the 
Casimir energy (solid, in red) and the total energy (dotted, 
in black) as function of the bag size. 



Ecas [9) 



1 



9 < tt/2 



— sm 



nag V4^ V I i^-O? 9>n/2 

+ C2sm^{9) + Cism'^{9) (22) 

+ ^6 sin6(0) + C8 sin«(0)), 



where the coefficients are: 

C2 = -0.13381, 
Cg = -0.01247, 



0.05085 
0.01241. 



(23) 



A few comments are in order regarding Ecas{9). First, 
Eq. (22) only contains the 9 dependent pieces of the 



Casimir energy; ^-independent terms, still proportional 
to do exist (see Ref. |21j) but do not effect our in- 
stanton calculation. Second, there is a sharp transition 
in Ecas at 9 — 7r/2, which can be traced to the energy 
of the lowest lying eigenmode crossing zero. Varying 9 
from < 7r/2 to > 7r/2, a bag eigenmode is dragged out of 
the Dirac sea and becomes a valence mode. This state 
is filled but its energy (valence quarks) is not contained 
in the Casimir energy, so it must be added separately. 
Once included, the valence modes make the interior bag 
energy completely smooth in 9. Likewise, the valence 
mode must be added separately to the baryon number 
of the Dirac sea (see Eq. (16)), keeping the net baryon 
number constant. 

Having reviewed the role of the bag action for static 
boundary conditions, we now need to see what hap- 
pens when the boundary conditions change as a function 
of time. In our 'bounce' solution describing Skyrmion 
unraveling, the boundary F{rhag) changes as a func- 
tion of time, while, simultaneously, the bag dynamics 
will influence the optimal path of this 'bounce' solution. 
As a first approximation, we assume the bag dynam- 
ics change adiabatically and simply change Ecas{9) 
Ecas{F{ri,agX{t))). In terms of the instanton calculation, 
this amount to changing the potential in Eq. [18} 



ViX)-^V{X)+EcasiFinagX)). 



(24) 



The new potential is shown below in Fig. [5] -the Casimir 
contribution is clearly important. We also point out that 
Ecas depends on r^ag in a relatively simple way, Ecas ~ 
''hag- Using the Casimir-improved V{X) we calculate the 
tunnehng exponential to be 6.7 x 10^^, 1.1 x 10~^, 1.5 x 
10-3 for rbag/R = 0.05,0.1,0.3. Adding the Casimir has 
shrunk the tunneling by roughly an order of magnitude. 

To see whether all consequences of time-dependent 
boundary conditions are captured by Eq 



(24) 



thorough investigation is necessary. As the detailed cal- 
culation provided in the next section shows, Eq (24) is 



not complete. An adiabatic calculation involves keeping 
the bag boundary static except for a short duration in 
which the bounce happens. Within the bounce, there are 
two special times which will prove to be important and 
whose effect is not included in Eq (24): at some time ti 



6 



I 



15 h 
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and we define {'n{T)\ to be the eigenstates of H{t). We 
have 



(n(T)l^l^(r)) + i?„(r)(n(r)|^(r)) = n (n(r)|V.(r)) 



dr 



;n(r)|V(r))+i?„(r)(n(r)|V'(T)>- 



(26) 
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FIG. 5: The energy profile including the Casimir energy 
(solid). The Skyrmion contribution, as shown in Fig. [s] is 
indicated by the dashed line. 



in the course of the unwinding, F{rhag) will become less 
than 7r/2, signifying that a mode from the Dirac sea has 
been lifted out. Similarly, as the Skymion rewinds, this 
same mode will dive back into the sea at some later time 
t2. This criss-crossing of F{rbag) = ""/S indicates that 
the system has a zero-mode, as shown in the cartoon of 
the Dirac sea in Fig. [2l This zero-mode has important, 
subtle implications^. 



(n(T)l^(T)) = K(n(r)|^(T)) 
Defining c„(t) = (n(r)|'0(T)), we obtain 

c„(r) -I- En{T)cn{T) - ^ (n(r) |TO(r)) c,„ (r) = KC„(r). 

m 

(27) 

If the boundary conditions arc changing slowly, the third 
term on the left-hand side is small and can be treated 
as a perturbation. We can rewrite the fermionic path 
integral as 



cexp 



d-r {4i (t) Aim (t) c„i (t) - 

cUr)K™(T)c„(r))l, (28) 



where Dnmir) = {dr + £^«(r))(5„m and Vnrn{r) = 
(n(r)|m(T)). Treating the first term in the exponent as 
the propagator and the second as a perturbation, the re- 
sult is 



A second consequence of time-dependent boundary det = det Dnm exp connected diagrams (29) 



condition absent in Eq (24 1 is that the Casimir energy 



depends on both F(ri,ag) as well as on its time derivative, 
F{fbag)- When manipulated into the bounce action, the 
F terms becomes A terms, and the coefficient of the 
term will play the role of a 'mass' for A. As such, it will 
affect the Skyrmion decay rate in a similar fashion to the 



The determinant of Dnm is easily evaluated because it 
is a disconnected set of one dimensional equations. The 
eigenfunctions c„(t) are 



K{X) term in Eq. (20 1 



c„(r) = exp 



j\ dT'En{T') 



(30) 



To determine k, we impose anti-periodic temporal 
boundary conditions, c„(r/2) + c„(— r/2) = 0: 



A. The non-static case 

For time-dependent 6{t) = F{ri,agX{T)), the result of 
the fermionic path integral is det(9T- -I- H{t)), where we 
have recast the time-dependent boundary conditions as 
a time-dependent Hamiltonian. In order to calculate the 
determinant we have to solve for the eigenvalues. Sup- 
pose 



kT-J"^ dr' Enir') ^ 2Tri(rn+^^ 
.27r(m+i) 



(31) 



dr' Enir') = iLUm + En 



where En = 



T/2 



T/2 



drEnir) (32) 



(9, + i7(r))|V'(T)) =A^|V(r)), 



(25) 



^ By the Atiyah-Patodi-Singer theorem |30H32j . if the system at 
any given instant has a zero mode solution, the full, time- 
dependent system will also exhibit a zero mode solution (see 
[33). 



In appendix |X] we provide an explicit calculation of 
det Dnm, determined by the product over all k. The re- 
sult, for time-dependent boundary conditions, is 

dctAim =cxp[-r^„,,] (33) 

n n 

Therefore, working to lowest order in an adiabatic ap- 
proximation, we see that the functional determinant is 
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the same as in the static case with the Casimir energy 
promoted to a function of time. Corrections to this result 
will be determined shortly, however, even at lowest order 
there is a subtlety incorporating Eq. ( 33 1 into the action 
for F{rbag)- For all n such that Enlj) doesn't change 
sign we have 



j dTEn{T)\ = JdT\EniT)\ 



(34) 



and there is no issue. However, as we noted regard- 
ing Eq. ( |22| ), the sign of the lowest energy eigenvalue 
changes as we cross F{rbag,T) ^ 9 = tt/2. In our case 
(the bounce) , Eq goes from negative to positive at ti and 
back at t2, leading to 



dr^o(T) 



dT\Eo{T)\-2 I rfr|So(T)|. (35) 



Therefore the function determinant becomes 

detD„„ = cxp { - Nc (^j drEcasir) + <^'' I^o(t)|^ }, 



The extra piece in Eq. ( 36 1 is exactly the valence term 



mentioned earlier. As we unwind the Skyrmion, we lift 
a negative state out of the sea and create the valence 
quarks. This means that we have yet another extra po- 
tential term in our instanton barrier. 



v{e)^E,asio) + 



Eo{e) 9 < tt/2 
9 > tt/2 



(37) 



In Fig.|6]the effect of this on the potential for k is shown. 
With the valence contribution, tunneling can no longer 
occur. The valence term makes the barrier in A impass- 
able, and the chiral bag is a stable configuration. It ap- 
pears Skyrmion unwinding is not possible, however the 
existence of the valence piece is intimately tied to the 
zero mode solution as we will show shortly 

Before discussing the zero mode, we calculate the per- 
turbative corrections to Eq. ( 36 ) . For this we need the 
propagator -D,7^ . Focusing on one energy level, we define 
two functions 



IpnfliT) = exp 

V'n,o(T) = exp 



dT'E,,{T') 

dr'E^ir') 



(38) 



These functions satisfy Dnm'ipn.o = 4'n,aDnm ~ 0. If 
£'„(±oo) < then we can use V'n.OiV'n.o to construct the 
states 



V'«,t'(t) = 'lpn,a{T)0{T' - t) 

4'n,r'{r)=i^^^oir)9{T-r'). 



(39) 



50 r 
40 ; 

^ 30- 
^ 20- 
10- 



0- 



10 



20 30 
X 



40 50 



FIG. 6: The energy profile including the Casimir and the 
valence quarks. The Casimir (dashed) and Skyrmion (dotted) 
contributions are the same as in Fig. |5] 



(40) 



These are normalizable functions for which it holds that 

Dnmlpnyir) = -'lpn,a{T')S{T - t') 
i^n,T'iT)Dnrn = ~^n,QiT')S{T - t') 

The inverse is solved through the following trick 

DD ^D'4jn,T2 = Snm'ipQ{Tl)lpQ{T2){Ti\D 
'4>m.TiDlpn,T2 = -S„m'ipoiT2)iJo{T2)0iT2 - Ti) 



(41) 



From which it follows 



exp 



drE^ir) 



9{t2 - Ti). 



(42) 



Similarily for positive energy levels we have 

(Ti|-D~™|r2) = (5,„„ exp -/ drEnir) 9{ti-T2). 

(43) 

Now we can systematically include the corrections by cal- 
culating Feynman diagrams order by order in insertions 
of Vnm- Because the vertex is of degree 2, the only con- 
nected Feynman diagrams are simple loops with n ver- 
tices, as shown in Figure. [T] The loop with one vertex 




+ 



FIG. 7: The first few perturbations of the fermion functional 
determinant. 

is zero, as we can always use rephasing freedom to set 
{n\n) — 0. The second order correction is given by 

- ^ J drdr' ^ Km(T)y„i„(T')Aim(T', T)i:)„i„(T, r'). 
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Due to the presence of 6'-functions in the propagators, the 
diagram is nonzero only if the energy levels n and m have 
opposite signs. Combining propagators (for i?„ > 0), we 
get 

(44) 

This propagator falls off exponentially in |r — r'| because 
of the difference in energy between n and m. Therefore, 
to a good approximation, we can just set r' — > r. Ex- 
panding VmniT') about r and keeping only the leading 
piece, the energy correction becomes 



dT 



\E,,{T)\ + \E^iT) 



(45) 



Using the chain rule, we have n = n'9. The energy cor- 
rection above therefore adds to the coefficient of the 9'^ 



term which, in analogy with Eq. ( 18 1, corresponds to the 
'mass' term of the instanton. 

Is this term ever important? Direct evaluation of 



Eq. ( 45 1 is arduous and beyond the scope of this work. 



however we can get some intuition from dimensional anal- 
ysis. From Eq. (45) we see that the non-adiabatic correc- 



tion is proportional to r^ag, while V^(A) is proportional 
to r'^^g. For bigger bag sizes, the potential barrier be- 
comes smaller, but the smaller potential is compensated 
by a bigger mass term. The effects of the non-adiabatic 
piece are exacerbated because the Skyrme contribution 
to K{\) (~ A~'^) becomes vanishingly small at large A 
- so any A-independent (or mildly A-dependent) piece 
of Eq. ( 45 ) will quickly become the dominant term. If 



K{\) goes to a constant, ^(A) alone will determine the 
tunneling rate and the suppression will get bigger. How- 
ever, without a full calculation of Eq. ( 45 1 we cannot be 



more quantitative. One option is to forget about the non- 
adiabatic corrections to K{X) and only use the Skyrmion 
contribution. This simplification limits us to small bag 
sizes, and really only gives us a lower limit on the size of 
the exponent. As a second option, we can parameterize 
the non-adiabatic contribution, 



K{\) K{^) + const X r^ag 



(46) 



To give a rough estimate, if the constant above is simply 
one, the exponential suppression becomes 



10- 



-35 



nag/R = 0.05 



2 X 10^16 ^^^^/^ ^ 0.1 
2 X 10-5 Tbag/R = 0.3 



However, before these numbers can be relevant, we need 
to somehow remove the stability exhibited in Fig. |6] 



B. The zero mode 



From Eq. (37) and Figure [6] it is clear that the Eq 
term from the valence quarks is what provides stability 



to the chiral bag system. The larger J^*^ dr|i?o(T)| is, 
the smaller the functional determinant and the larger the 
barrier inhibiting Skyrmion decay. How can Skyrmion 
decay occur in the face of this barrier? The answer lies 
in recalling that zero modes (or approximate zero modes) 
also lead to suppressed functional determinants. 

In a system with zero modes, despite the fact that 
the functional determinant is vanishing, it is well known 
that correlations functions of operators can be nonzero. 
This fact is most easily proven by adding sources to the 
fermion Lagrangian, as functional derivatives with re- 
spect to the sources cancel the zeros in the determinant. 
To account for these non-vanishing correlation functions, 
a semi-local operator - the 't Hooft vertex [28l[29] - must 
be added to the action. 

Proceeding by analogy, when dT\EQ{T)\ is large 
there are approximate left and right zero modes of D. 
The existence of the zero modes makes the effective 
Skyrmion action zero unless sources are added. Including 
sources, the result of the fermionic path-integral Z[J, 
is 

exp -Nc (^j drE^Ur) + ^ ' dT\Eo{T)\^ + f D^\j 

(48) 

From Eq. ( 42 ) , we can see that functional derivatives with 
respect to the sources are accompanied by D^^ ^ which 
has exactly the right exponent to cancel the (otherwise 
large) j^^ dT\Ea{T)\ piece. Stated another way, in order 
for the Skyrmion to unwind, there must be some (micro- 
scopic) source for baryon- number violation present, such 
as the higher-dimensional operator 

-■4)^11) j..^NcX = Ae'^-- '^ipi^ljj...ipNc- 

(49) 

Here x represents some lepton and M is the scale sup- 
pressing these operators; for convenience, we combine ev- 
erything into a (fermionic) coefficient A. The need for 
these operators is no surprise - we wouldn't be consid- 
ering proton decay at all if they were forbidden! To see 
how the operator in Eq. (49 1 generates the Dq^ needed 
to cancel the Eq term, rewrite the interacting theory as 



exp • 



{-/(*(^)"-^*(^)")H"''- <»> 



Expanding out Z[J,J^, the only nonzero term has Nc 
of both J and J^, and therefore comes with Nc copies 
(47) of Dq^ - just the right amount to cancel the Eq term in 



the functional determinant. The relevant term being 



Z[J,f] 3 exp 



-Nc / dtEcasit) 



{jHh)J{t,)fo, 



(51) 

The J'^ destroys a particle at ti - when its energy be- 
comes positive - and J creates it when its energy becomes 
negative. The fermionic sources are cancelled when we 



apply the functional derivatives (in Eq. (50)) correspond- 



ing the baryon number violating interaction. 
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The resulting picture of proton decay becomes thus of 
a meta-stable chiral bag with only the negative states 
filled. There is an instanton bounce solution, where the 
Skyrmion unwinds pulling the valence quarks out of the 
vacuum at which point they are destroyed by the baryon 
violating operator. The bounce solution can be calcu- 
lated by adding the Casimir energy to the potential, as 
postulated in the previous section. The bounce comes 
with the prefactor A^^ A. Therefore, the width of the pro- 
ton is given by 



(52) 



The AA''' gives the precise perturbative suppression of 
Eq. [l]as it should. We find, however, the extra exponen- 
tial suppression of the hadronic matrix element coming 
from the unwinding. In the above, we have ruthlessly 
suppressed any other factor. The calculation as it stands 
will not give a precision result, however the presence of 
the exponential suppression in the hadronic matrix ele- 
ment is our main result, and it can be very relevant. 



VI. DISCUSSION 

In this work we have provided a rigorous instanton 
calculation for the decay of a "baryon" in the context of 
the chiral bag model. The chiral bag model presents a 
model of the baryon where there are no valence quarks. It 
is a model where the valence quarks attract anti-quarks 
from the surrounding Dirac sea like a charge in a di- 
electric. In doing this it generates a topological twist 
in the chiral phase of the condensate, described by the 
Skyrmion of the pion field. The binding of the valence 
quarks with the anti-quarks lowers their energy, driving 
them into the Dirac sea of negative states. If this is 
an accurate description of the proton a non-perturbative 
suppression is to be expected. Removal of the valence 
quarks leaves a state with higher energy and therefore it 
must proceed through a tunneling process. 

The calculation we have done remains sensitive to the 
bag radius, and therefore does not obey the 'Cheshire Cat 
principle' as nicely as one may have hoped. The chiral 
bag is, however, just a leading order factorization model 
of hadronic QCD: the non- interacting free quark theory 
within the bag cannot possibly generate the spontaneous 
symmetry breaking necessary for the pion field, so the 
bag size cannot be arbitrarily large. Similarly, the pion 
field description also breaks down for small bag size, due 
to the running of the coupling. Given the limitations 
of the chiral bag, dependence on rbag it not surprising. 
Numbers aside, our main qualitative point is that this 
model gives the right picture of what a proton is, and 
that therefore a tunneling suppression is to be expected. 
If a more sophisticated calculation including NLO effects 
results in a regime in which the calculation stabilizes, this 
would certainly strengthen our result. 

Applying the suppression we find to simple GUT mod- 
els has profound implications. An additional ^ 10~^ 



suppression would mean the unification scale could be 
lowered by an order of magnitude for dimension-six 
baryon-number violation, or two orders of magnitude 
for dimension-five baryon number violation. Further- 
more, 10^^ should be viewed as an upper limit. In- 
cluding a parametrization of the non-adiabatic terms (see 
Eq. (47)), we found 10~^^ suppression was perfectly rea- 
sonable, implying a three (six) order of magnitude drop 
for dimension six (five) baryon number violation. A more 
exact value for the exponent requires a rigorous calcula- 
tion of Eq. (45). A change in the required GUT scale 



of this order would certainly resuscitate several scenar- 
ios [SiMO] ! According to our calculation other tests of 
baryon number violation, such as neutron-antineutron 
oscillation |1T| - H5] . should also be highly suppressed. 

Throughout our calculations we have completely ig- 
nored any mass for the interior quarks. Adding in a mass, 
the energy levels shift, and if the shift is large enough the 
mass can effectively prevent valence quarks from diving 
in the vacuum. If the valence quarks retain positive en- 
ergy, they can decay immediately and do not need to be 
lifted by unwinding the pion field. The baryon is then 
an ordinary bound state. The role quark masses play in 
this calculation may also explain the apparent disagree- 
ment between our result and estimations of proton decay 
matrix elements based on lattice QCD [J'-l]. The sup- 
pression in the chiral bag model comes primarily from 
the chiral symmetry and its twisting. Exact chiral sym- 
metry is a difficult regime to probe on the lattice due 
to difficulties with chiral symmetry of fermions and with 
fitting the Compton wavelength of low mass states into 
the finite lattice volume. Lattice results therefore depend 
on extrapolation of quark mass and system volume into 
the physical regime. If the dependence of proton decay 
matrix elements on lattice artifacts is different than for 
the more conventional observables (meson masses, etc.), 
extrapolations, motivated by conventional observables, 
would be inappropriate and may explain the apparent 
differences between our result and the lattice. The mass- 
dependence in our result is a subtle issue and deserving 
of additional study. 



Appendix A: Deriving the Casimir energy 



The eigenvalue spectrum is given by 



27r (m + i 
f 




(Al) 



This spectrum is symmetrical around the real axis, so we 
can pair every conjugate pair and we have 



Using zeta-function regularization we have 
1 d 



det D = exp 



2 ds 



(A2) 



(A3) 
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In the limit T — ^ oo we can replace — >■ ^ J duj, we 
obtain 



1 



)/r(s) (A4) 



Because of the T[s] = r[s + l]/s we see that the above 
expression is at s = 0. Therefor the only way to get 
a non-zero values is if the differential operator is applied 
to kill s. We are thus left with 



detL' = exp[-J^|^„ 



This can be written as 

deti:i = exp[-£:casT] 
where the Casimir energy Ecus is 



ZTT Z 



(A5) 
(A6) 



s=0 



(A7) 



- exactly the zeta-function regularization of the vacuum 
energy. 



Acknowledgments 



We thank Paddy Fox, Roni Harnik, Chris Hill, Andreas 
Kronfeld and Ethan Neil for many valuable conversations 
and acknowledge Joachim Kopp for collaboration during 
the early stages of this work. GS acknowledge Piet Mul- 
ders for a useful discussion. AM and GS are supported 
by Fermilab operated by Fermi Research Alliance, LLC 
under contract number DE-AC02-07CH11359 with the 
US Department of Energy. 



H. Georgi, S. L. Glashow, Phys. Rev. Lett. 32, 438-441 

(1974) . 

N. Tsutsui et al. [CP-PACS Collaboration and JLQCD 
Collaborations], Phys Rev. D 70, 111501 (2004) 
arXiv:hep-lat /0402026" . 

Y. Aoki et al. [RBC-UKQCD Collaboration], Phys. Rev. 
D 78, 054505 (2008) arXiv:0806.1031 [hep-lat]]. 
P. Cooney, http://www.era.lib.ed.ac.uk/handle/1842/4042 
M. Claudson, M. B. WiseTLTXWlT, NucirPFysTBTMr 
297 (1982). 

R. Jackiw, C. Rebbi, Phys. Rev. D13, 3398-3409 (1976). 
V. Vento, M. Rho, E. M. Nyman, J. H. Jun, G. E. Brown, 
Nucl. Phys. A345, 413 (1980). 

J. Goldstone, F. Wilczek, Phys. Rev. Lett. 47, 986-989 
(1981). 

J. Goldstone, R. L. Jaffe, Phys. Rev. Lett. 51, 1518 
(1983). 

G. E. Brown, A. D. Jackson, M. Rho, V. Vento, Phys. 

Lett. B140, 285-289 (1984). 

P. J. Mulders, Phys. Rev. D30, 1073 (1984). 

L. Vepstas, A. D. Jackson, A. S. Goldhaber, Phys. Lett. 

B140, 280-284 (1984). 

L. Vepstas, A. D. Jackson, Phys. Rept. 187, 109-143 
(1990). 

T. H. R. Skyrme, Proc. Roy Soc. Loud. A260, 127-138 
(1961). 

T. H. R. Skyrme, Nucl. Phys. 31, 556-569 (1962). 

A. Chodos, R. L. Jaffe, K. Johnson, C. B. Thorn, 

V. F. Weisskopf, Phys. Rev. D9, 3471-3495 (1974). 

T. A. DeGrand, R. L. Jaffe, K. Johnson, J. E. Kiskis, 

Phys. Rev. D12, 2060 (1975). 

T. Inoue, T. Maskawa, Prog. Theor. Phys. 54, 1833 

(1975) . 

C. G. Callan, Jr., R. F. Dashen, D. J. Gross, Phys. Rev. 
D19, 1826 (1979). 

G. E. Brown, M. Rho, Phys. Lett. B82, 177-180 (1979). 
G. E. Brown, M. Rho, V. Vento, Phys. Lett. B84, 383 



(1979). 

[22] S. Nadkarni, H. B. Nielsen, I. Zahed, Nucl. Phys. B253, 
308 (1985). 



[23 
[24 

[25 
[26 
[27 

[28 
[29' 
[30 

[31 

[32 

[33 
[34 

[35 

[36 
[37 

[38; 

[39 

[40 
[41 

[42; 



M. Rho, Phys. Rept. 240, 1 (1994). [hep-ph/ 0206003 . 

M. De Francia, H. Falomir, E. M. Santangelo, Phys. Lett. 

B371, 285-292 (1996). hep-ph/9507347 . 

E. Witten, Nucl. Phys. B223, 433-444 (1983). 

E. Witten, Nucl. Phys. B223, 422-432 (1983). 

G.H. Derrick, J. Mathematical Phys. 5 (1964), pp. 

12521254 

G. 't Hooft, Phys. Rev. D14, 3432-3450 (1976). 

G. 't Hooft, Phys. Rev. Lett. 37, 8-11 (1976). 

M. F. Atiyah, V. K. Patodi, I. M. Singer, Math. Proc. 
Cambridge Phil. Soc. 77, 43 (1975). 

M. F. Atiyah, V. K. Patodi, I. M. Singer, Math. Proc. 

Cambridge Phil. Soc. 78, 405 (1976). 

M. F. Atiyah, V. K. Patodi, I. M. Singer, Math. Proc. 

Cambridge Phil. Soc. 79, 71 (1980). 

M. Nakahara, Geometry, topology and physics (2003) 

J. R. Ellis, D. V. Nanopoulos, S. Rudaz, Nucl. Phys. 

B202, 43 (1982). 

P. Nath, A. H. Chamseddine, R. L. Arnowitt, Phys. Rev. 
D32, 2348-2358 (1985). 

P. Nath, R. L. Arnowitt, Phys. Rev. D38, 1479 (1988). 
J. Hisano, H. Murayama, T. Yanagida, Nucl. Phys. 
B402, 46-84 (1993). hep-ph/9207279 . 

H. Murayama, A. Pierce, Phys. Rev. D65, 055009 (2002). 
|hep-ph/0108104 . 

T. Fukuyama, A. Ilakovac, T. Kikuchi, S. Meljanac, 
N. Okada, Eur. Phys. J. C42, 191-203 (2005). [hep^ 
ph/0401213. 

B. Dutta, Y. Mimura, R. N. Mohapatra, Phys. Rev. Lett. 

94, 091804 (2005). hep-ph/0412105 . 

V. A. Kuzmin, Pisma Zh. Eksp. Teor. Fiz. 12, 335-337 

(1970). 

R. N. Mohapatra, R. E. Marshak, Phys. Rev. Lett. 44, 
1316-1319 (1980). 



11 



[43] L. N. Chang, N. P. Chang, Phys. Lett. B92, 103-106 

(1980). 



